home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zgttrf.z / zgttrf
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))                                                          ZZZZGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZGTTRF - compute an LU factorization of a complex tridiagonal matrix A
  10.      using elimination with partial pivoting and row interchanges
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZGTTRF( N, DL, D, DU, DU2, IPIV, INFO )
  14.  
  15.          INTEGER        INFO, N
  16.  
  17.          INTEGER        IPIV( * )
  18.  
  19.          COMPLEX*16     D( * ), DL( * ), DU( * ), DU2( * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      ZGTTRF computes an LU factorization of a complex tridiagonal matrix A
  23.      using elimination with partial pivoting and row interchanges.
  24.  
  25.      The factorization has the form
  26.         A = L * U
  27.      where L is a product of permutation and unit lower bidiagonal matrices
  28.      and U is upper triangular with nonzeros in only the main diagonal and
  29.      first two superdiagonals.
  30.  
  31.  
  32. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      DL      (input/output) COMPLEX*16 array, dimension (N-1)
  37.              On entry, DL must contain the (n-1) subdiagonal elements of A.
  38.              On exit, DL is overwritten by the (n-1) multipliers that define
  39.              the matrix L from the LU factorization of A.
  40.  
  41.      D       (input/output) COMPLEX*16 array, dimension (N)
  42.              On entry, D must contain the diagonal elements of A.  On exit, D
  43.              is overwritten by the n diagonal elements of the upper triangular
  44.              matrix U from the LU factorization of A.
  45.  
  46.      DU      (input/output) COMPLEX*16 array, dimension (N-1)
  47.              On entry, DU must contain the (n-1) superdiagonal elements of A.
  48.              On exit, DU is overwritten by the (n-1) elements of the first
  49.              superdiagonal of U.
  50.  
  51.      DU2     (output) COMPLEX*16 array, dimension (N-2)
  52.              On exit, DU2 is overwritten by the (n-2) elements of the second
  53.              superdiagonal of U.
  54.  
  55.      IPIV    (output) INTEGER array, dimension (N)
  56.              The pivot indices; for 1 <= i <= n, row i of the matrix was
  57.              interchanged with row IPIV(i).  IPIV(i) will always be either i
  58.              or i+1; IPIV(i) = i indicates a row interchange was not required.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))                                                          ZZZZGGGGTTTTTTTTRRRRFFFF((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      INFO    (output) INTEGER
  75.              = 0:  successful exit
  76.              < 0:  if INFO = -i, the i-th argument had an illegal value
  77.              > 0:  if INFO = i, U(i,i) is exactly zero. The factorization has
  78.              been completed, but the factor U is exactly singular, and
  79.              division by zero will occur if it is used to solve a system of
  80.              equations.
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.